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0.  Introduction 


This  i>aper  provides  a stability  study  of  a vide  class  of  difference 

• ■ 

[ approximations  for  a hyperbolic  mixed  lnltlal>boundary  value  problem  in 

1 

\ the  qviarter  plane  x > 0,  t > 0.  The  approximated  differential  system 

is  U4.  s Au_  vhere  A is  diagonal,  and  the  inflow  and  outflow  unknowns 

w X 

interact  at  the  boundary.  For  the  difference  approximation  we  consider 


general  dissipative  explicit  two>level  schemes,  with  general  boundary 
conditions  which  determine  the  boundary  values  in  terms  of  outflow  values 
at  interior  points.  We  show  that  such  boundary  conditions  may  have 
arbitrary  degree  of  accuracy.  This  is  included  in  Section  1. 

In  Section  2 we  begin  to  discuss  stability.  We  show  that  the  entire 
approximation  is  stable  if  and  only  if  the  scalar  coinponents  of  its  out- 
flow part  are;  thiis  reducing  the  stability  question  to  that  of  a scalar 
outflow  problem.  From  that  point  on,  our  aim  will  be  to  obtain  easily 
checkable  sufficient  stability  criteria  for  the  reduced  problem.  All  our 
results  sore  scheme-independent  and  are  given  exclusively  in  terms  of  the 
outflow  boundary  conditions.  The  only  such  result  that  we  -know  of,  is 
due  to  Kreiss  ([3]j  see  also  [1])  who  proved  that  for  dissipative 
schemes,  boundary  extrapolation  always  maintains  stability. 

In  the  remainder  of  Section  2 we  state  our  main  result  and  consider 
seversl.  examples.  The  main  result  is  for  the  case  where  the  outflow 
boundary  conditions  sa-e  translatory,  i.e.,  determined  at  all  boundary 


points  by  the  same  procedure.  The  result  states  that  if  the  outflow 


boundary  conditions  are  generated  by  a solvable  stable  scheme,  then  the 
entire  approximation  is  stable,  independently  of  the  interior  scheme. 

The  examples  considered  show  that  if  the  outflow  boundary  conditions  are 
generated  by  oblique  extrapolation,  by  explicit  or  implicit  Euler  schemes. 


2 

or  by  the  Box-Scheme,  then  overall  stability  Is  assured. 

In  Section  3 we  derive  e preliminary  scheme-independent  result, 
which  Is  an  e]q>llclt  Interpretation  of  Xrelss'  determlnantal  stability 
criteria  [i)'].  Ihis  result,  lAich  seems  to  be  of  independent  Interest, 

Is  used  In  Section  4 to  derive  our  stability  criteria  for  the  cue  of 
traulatoxy  boundary  conditions. 

The  theoretical  buls  for  our  work  was  given  by  IQrelss  [4]  and  by 
Oustafsson,  Krelss  and  Sundstrbm  [2].  We  usume  that  the  reader  Is 
familiar  with  these  papers. 

In  a forthcoming  paper  we  discuss  the  extension  of  our  results  to 
nondlsslpatlve  and  multi-level  schemes. 

We  are  grateful  to  BJom  Engqulst  for  many  helpful  discussions. 


1.  The  Difference  Apiaroxlmatlon 

Consider  a first  order  hyperbolic  system  of  partial  differential 
equations 


(1.1a) 


du/dt  > A2u/dx,  X > 0,  t'  > 0, 


idiere  u ■ (u^^^(x,t),...,u^**\x,t))'  is  the  transposed  vector  of  un- 
knowns, and  A is  a constant  nxn  Hermltlan  matrix  of  the  form 

A 


, A^  < 0,  A^  > 0, 


0 A 

Without  restriction  we  may  usume  that  A is  diagonal. 


Hie  solution  of  (1.1a)  Is  uniquely  determined  If  Initial  values 


(1.1b) 


u(x,0)  - f(x),  x>0. 


Onmr  ■ **  w ■■ 


■>■>1  'ir'C* 


DO 
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and  boundary  conditions 


(1.1c) 


U^(0,t)  - Su^^(0,t)  + g(t),  t > 0, 


are  prescribed.  Here  the  partition 

* 

uI.(uW uW)' 

corresponds  to  that  of  A,  and  S is  an  i x (n  - /)  rectangular  matrix. 

To  solve  the  initial-boundary  value  problem  (1)  by  a difference 
approximation  we  introduce  a mesh-size  h*Ax>0,  hEAt>0,  such  that 
X B k/h  = constant.  Using  the  standard  notation  v^(t)  = v(vh,t),  we 
approximate  (1.1a)  by  a consistent  two-sided  difference  scheme  of  the  form 


(1.2a) 


\(t  + k)  = Qv^(t),  V = r,r+l,..., 


Q = S » 


d=-r 


AjE-, 


Ev  = V 
V v+1* 


with  initial  values 


(1.2b) 


V (0)  - f , 
V * v' 


V c 0^192^...  . 


The  Aj  are  fixed  diagoneil  nxn  matrices  depending  on  A and  on  X 
sack  that  A , A are  nonsingular. 

-r»  p 

Ihrou^out  the  paper  we  assume  that  scheme  (1.2a)  is  dissipative,  i.e., 
for  some  6 > 0 and  positive  Integer  the  eigenvalues  6(|)  of  the 
amplification  matrix 


0,(0  = h A.e^^^ 

d— r 


-IT  < I < TT, 


|e(0l<i-BUl^“ 


satisfy 


! 


Since  the  Aj  are  diagonal,  dlssipativity  guarantees  that  the  scheme 
(1.2a)  is  stable. 

In  order  to  uniquely  deteimlne  the  solution  of  (1.2),  we  must  spec- 
ify, at  each  time  step,  the  r boundary  values  'v^(t),  4 « 0,1,..., r-1. 
For  the  outflow  unknowns  we  do  It  by  boundary  conditions  of  the  form 


(1.3) 


Z)  (t  + k)  - E (t). 

mJ  4+J'  J.0  4J  4+j' 


being  fixed  diagonal  (n  - je)  x (n  - i)  matrices.  For  the  Inflow  part 
we  use  the  physiced  boundary  condition 


(l.l»a) 


vj(t)  - SvJ^(t)  + g(t) 


together  with  r«l  additional  conditions  of  the  form 


(1.4b) 


4 “ 1, . . . , r-1. 


where  are  fixed  i x (n  - f)  matrices  and  the  g^(t)  depend  on  h 

and  on  g(t). 

It  Is  well  known  that  using  conditions  of  the  general  form  (l.3)> 
one  can  achieve  at  the  boundary  arbitrary  degrees  of  accuracy.  We  note 
that  this  Is  true  also  for  conditions  of  type  (1.4b).  In  fact.  If  accuracy 
of  order  d Is  desired,  then  using  the  differential,  system  and  (1.1c), 
we  find  that  a Taylor  ercpanslon  of  a smooth  solution  of  (1.1)  yields 


ul(t)  - s ui(o.t)  * 

. h {aI)-J  4 (o.t)  ^ 

j-0  •'*  afJ 

. £ (aI)-^  is  4 u“  (o.t)  . 4 " 'X"**") 

J-0  ''*  at*'  dt*' 

. S (A^)-J  [S(a”)^  4 («.*)  - 4 «(*)>  * 

j-o  ax^  dt^ 

Thus, (1.4b)  follows  upon  approxixnating  a^/ax^u^^(0,t)  by  linear  combina- 
tions of  u^^  ('4)>.. . ,u^^(t)  of  the  ri^t  accuracy. 

For  exan^le,  if  v^(t)  is  required  to  second  order  of  accuracy,  we 


may  use 


3uj"(t) . kup(t)  - >4^(t)  jv^fo « 


II/.N  . ..II/ 


to  obtain  a second  order  accurate  boundary  condition 


vj(t)  - DQvJ^(t)  + D3^vJ^(t)  + D2V^^(t)  + g^^Ct), 


where 


®0  “ 4) 'I  “l'*‘¥  Mg,  Dg  - + J Mg, 


Mj  - (A^)’J  S(A^^)**, 


i •=  0,1,2, 


g3^(t)  - g(t)  + h(A^)‘^  g(t)  + ~(A^)’^  g(t). 

2.  Statement  of  Main  Result  and  Examples 

The  difference  approximation  is  completely  defined  now  by  the  dissi 
pative  (stable)  scheme  (1.2)  together  with  the  boundary  conditions  (1.3) 
(1.4),  and  we  raise  the  question  of  overall  stability  in  the  sense  of 
Oustafsson  et  al.  ([2],  Definition  3.3)* 


i 

J 

■ 

I 


6 


\ 

i 

t 

\ 


I 


f 


Since  the  are  diagonal,  ee  can  split  the  scheme  (1.2a)  into 
its  inflow  ai^  outflow  parts, 

(2.1)  ^(^  ^ A?  ''^^('^)»  ^ “ r,r+l, ..., 

^ d—r  J . 

and 

(2.2)  (t  + k)  - h a“  vJJjCt),  . V - r,rfl,,,., 

J*“r 

idiere 

r < d < P, 

corresponds  to  the  partition  of  A.  We  immediately  see  that  the  outflow 
problem  (2.2),  (1.3)  is  self  contained,  while  the  inflow  problem  (2.1), 
(l.U)  depends  on  the  outflow  part  only  to  the  extent  that  the  outflow 
ccapatatlons  provide  the  inhomogeneous  boundary  values  in  (1.4).  Thus, 
stability  of  the  entire  approximation  is  equivalent  to  the  following 
separate  questions: 

(a)  Stability  of  the  inflow  problem  (2.1)  with  inhomogeneous  boundary 

values. 

(b)  Stability  of  the  outflow  problem  (2.2)  (1.3). 

Since  the  stability  definition  (3.3)  of  [2l  gives  bounds  for  the 
inhomogeneous  boundary  values,  it  suffices,  for  the  outflow  problem,  to 
consider  homogeneous  boundary  values.  The  Aj  are  diagonal,  hence  the 
problem  splits  into  i independent  dissipative  approximations  with 
homogeneous  boundary  values.  It  was  shown  by  Krelss  ([3],  Theorem  3)  that 
such  scalar  approximations  are  stable  Independently  of  the  basic  scheme. 
Thus  the  Inflow  problem  is  stable  and  it  remains  to  consider  (b).  Since 
the  Aj^  of  (2.2)  and  the  of  (1.3)  are  diagonal,  the  outflow  prob- 

lem splits  as  well. 
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it  the  n - 1 scalar  components  of  Its  outflow  part  are  stable.  So, 
frcB  now  on,  we  may  restrict  the  stability  discussion  to  the  following 
scalar  case  where  we  approximate  an  inltial-V|^ue  problem 


(2.3) 


- a |H,  a > 0,  X > 0,  t > 0,  u(x.O)  - f(x), 


by  a consistent  dissipative  scheme 


(2.4a) 


\(t  + k)  - Qv^(t), 


Q-  £ a.E^, 
J-r  ^ 


together  with  initial  values 


V - r,r+l,..,. 


• Vl- 


(2.lrt>) 


V (O)  « f , 
V ' V 


and  boundary  conditions 


V ■ 0,1,2,..., 


d5)  " 0,...,r.l. 

Here  the  a^  and  the  c^  are  constants  depending  on  a and  on  1. 

Our  purpose  is  to  provide  easily  checkable  sufficient  stability 
criteria  fbr  the  approximation  (2. 4)  (2. 5)  idilch  depend  entirely  on  the 
boundary  conditions  (2.5).  Particularly  effective  criteria  of  that 
nature  are  obtained  when  the  boundary  conditions  are  translatory.  i.e. , 


of  the  foiv 


M - 0,...,r-l, 


idiere  the  c^  are  Independent  of  u.  In  order  to  state  our  main  result 
we  Introduce  the  boundary-scheme 
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Y^(t  + k)  - v^(t),  V - 0,+l,+2,..., 

(2.7)  , 

t(®)  ■ T^®^(E)  - E or  - 0;i,  Ev  . ▼ 

J-0  ^ V v+l 

idilcb  generates  the  boundary  conditions  (2.6)  upon  restriction  of  v to 
the  values 

We  call  the  boundary  scheme  stable  if  it  is  stable  vhen  applied  to  all 
grid  points  vh,  - • < v < «.  This  leads  to 

THEOREM  2.1  (The  Main  Theorem).  Approximation  (2.U)(2.6)  is 
stable  if  the  boundary-scheme  (2.7)  is  stable  and  if 

(2.8)  T^^^k)  e E c$^)  / 0 VUI  < 1. 

J-0  ^ 

The  stability  criterion  in  this  theorem  is  independent  of  the  basic 
scheme.  The  proof,  as  veil  as  other  scheme -independent  results  for  the 
translatory  case,  are  given  in  Section  U. 

Often,  the  boundary- scheme  is  knovn  in  advance  to  be  stable.  Thus, 
in  applying  Theorem  2.).,  it  only  remAlns  to  verify  the  solvability- 
condition  (2.8). 

If  the  boundary  conditions  are  explicit,  i.e. , of  the  form 

(2.9)  T^(t  + k)  - Cj  v^^j(t),  n - o,...,r-i, 

then  T^^^(k)  c 1.  Hence  (2.8)  is  automatically  fulfilled  and  Theorem  2.1 
reduces  to 

COROUARY  2.2.  If  the  boundary  conditions  (2.9)  are  generated  by  a 
stable  scheme,  then  approximation  (2.4)(2.9)  is  stable. 

EXAMFIE  1.  Detenalne  the  boundary  conditions  by  oblique  extrapolation 





(2.W) 


V*  * ‘ 


^ • Oy • • • 


Clearly,  tBe  generating  boundary- scheme  Is  ejqpUcit  and  uncondltionedly 
stable;  so  by  Corollary  2.2  the  stability  of  (2.U)(2.10)  is  assured. 

note  that  (2.10)  is  inconsistent  vith  the  differential  equation 
(2.3)  unless  Xa  K 1 in  which  case  (2.10)  coincides  with  our  next  example. 

EXAMPI£  2.  Let  the  boundary  conditions  be  generated  by  the  rl^t- 
sided  explicit  Euler  scheme,  i.e.. 


(2.U) 


v^(t  + k)  - v^(t)  + \a(v^^^(t)  - v^(t)],  ix  - 0,...,r-l. 


Since  the  basic  scheme  (2.ba)  must  satisfy  the  Courant-Frledrlchs-Levl 
condition  Xa  < 1,  Euler's  scheme  is  stable  too,  and  by  Corollary  2.2 
the  stability  of  (2.1).)(2.11)  follows. 

EXAMPLE  3.  Take 

(2,12)  v^(t  + k)  - Xa[v^^j^(t  + k)  - v^(t  + k)  ] = v^(t) , 

0,,..,r— 1, 

idiich  is  generated  by  the  right-sided  unconditionally  stable  Isqpllcit 
Euler  scheme.  To  comply  with  Theorem  2.1,  we  must  have 

T^^)(k)  e 1 - Xa(K  - 1)  / 0 V|k1  <1. 

Since  Xa  > 0,  then  k with  |k|  < 1 gives 

Re  T^^^(k)  » 1 + Xa(l  - Re  k)  > 1, 
so  (2. 3) (2. 12)  is  stable. 


EXAMPIE  U.  We  use  the  unconditionally  stable  Box-Scheme  to  generate 


ID 


(2.13)  v^(t  + k)  + + k)  - Xa(v^^^(t  + k)  - v^(t  + k)] 

Since 

Re  T<")(«)  . 1 + Re  K + Xa(l  - Re  k)  > 0 V Itl  < 1, 

then  by  Theorem  2.1  again,  (2. h) (2.13)  1b  stable. 

m view  of  the  previous  stability  discussion,  Examples  1>^  iiqply  that 
if  the  boundary  conditions  (1.3)  8nre  generated  by  oblique  extrapolation, 
explicit  or  iniplicit  ri^t-slded  Euler  schemes,  or  by  the  Box-Scheme,  then 
the  entire  approximation  (1.2)  - (1.4)  is  stable. 

The  boundary  conditions  in  Exanqples  1-4  were  studied  by  Gustafsson 
et  al.  [2]  and  by  SkbUermo  [3]  in  combination  with  specific  3-Point  basic 
schemes. 

3.  A Preliminary  Determlnental  Criteria 

In  order  to  investigate  the  stability  of  (2.4),  (2.3)  we  introduce 
the  space  jt2(b)  ol*  &11  ffrid  functions 

^ “ KCo  ® ^ l\l^  < -• 

VbO 

We  write  the  approximation  in  operator  form 

v(t  + k)  «:  Ov(t), 

where  G:  je2(M  uniquely  determined  by  (2.4a)  and  (2.3). 

Following  Krelss  ([4],  Main  Theorem)  we  will  show  that  0 has  no  eigen- 
values z with  |z|  >1.  This  will  suffice  to  assure  stability  both  in 
the  sense  of  [4]  and  in  the  sense  of  Definition  3.3  of  [2]. 


n 


To  check  the  eigenvalues  of  G ve  must  adopt  Kreiss*  recipe  in 
[4]:  If  z ifith  |z|  > 1 is  an  eigenvalue,  then  for  some  nontrivial 
g .6  ^8  * zg.  Thus,  by  (2.4a)  and  (2.5),  g must  satisfy  the 

resolvent  equation 


(3.1) 


and  the  boundary  relations 


r,r+l,.. . , 


(3.2) 


Z S c(})g  . . = E c%  . 
jcO  ^ 4'*"d  J=0 


(i  s 0, . . . , r-1. 


The  most  general  solution  of  (3.1)  in  l2(^)  be  written  as 

k 

(3.3)  g^  = E E o^P«ft(v)KX,  V > 0 . 


o - “ a* 

Ofcl  p=0 

Here  the  •'^r  * 1 < ® < 4,  are  the  distinct  roots  of  the  character- 


istic equation 


(3.4) 


V 4 

E a.K**  - z = 0, 


which  satisfy  0 < 1k^|  <1,  each  with  multiplicity  m^  ® “q(z);  the 
P^(v)  are  arbitrary  polynomials  in  v with  deg[P^(v)]  = Pj  and  the 
axe  free  coefficients  yet  to  be  determined.  Since  (2.4a)  is  dissipat- 
ive and  consistent  with  (2.3),  we  \ise  Lemma  2 and  part  of  the  proof  of 
Lemma  7 of  [4]  to  find  that  for  z with  |z|  >1,  equation  (3.4)  has 
precisely  r roots  k with  0 < |ic|  <1.  Thus  g^  of  (3.3)  depends  on 
r pexameters.  Substituting  (3-3)  in  (3.2),  we  obtain  a linear  homo- 
geneous system  of  r equations  with  the  r unknowns  The  system 

Op 

is  of  the  form  Ja'  » 0 where  J * j(z)  is  an  rxr  matrix  and  a is 
the  unknown  vector.  This  yields 


12 


mWi  3.1  (Lemma  3,  [4]).  z vlth  |z|  > 1 is  an  eigenvalue 
of  0 If  and  only  If  det  j(z)  / 0. 

Qolng  through  the  above  process,  we  make  a particular  choice  of  the 
polynomials  that  leads  to  an  explicit  expression  of  J idiich 

later  proves  useful.  We  choose 


(3.5) 


Inserting  (3.5)  in  (3.3)  ®nd  then  in  (3.2)  we  obtain 


(3.6)  Z Z Z [zoW  - 


Obl  JbO 


U • 0,..., r™!, 


idiich  constitutes  the  system  Jo*  ^ 0. 

At  this  point  we  associate  with  the  boundary  conditions  (2.5)  a set 
of  polynomial  boundary-functions 


(3.7)  b^(z,K)  . tz=W  . 

Since 

imO  ^ ^ ^ 

then  system  (3.6)  may  be  written  as 


h ■ 0,%.., r-l. 


(3.8) 


k (z,kJ 


Oh.1  p-O  bK' 


p r*l* • 


Thus  the  coefficient  matrix  J takes  the  form 


(3.9a) 


J ■ [B(z, • ,B(z, K|^,mj^)  ] 


where  the  B(z,K^,m^),  1 < a < k,  are  rxm^  blocks  given  by 
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(3.9b) 


bo(t,K) 

1 

X 

1 

Vi 

bQ(*.0 

b3^(z,K) 

a 

a 

» ^ 

bi(z,K) 

a 

a 

. 0 



Sk 

bj^(»,K) 

a 

a 

a 

a 

b (*» '^) 

a 

by.i(Z,K) 

b -.1  (*»•') 

II 

Defining  a partition  of  r to  be  any  set  of  positive  integers 
^ich  satisfies  + . . . + bIj^  » r,  we  state 

THEOREM  3.2.  Approximation  (2.4)(2.5)  Is  stable  if  for  every  z 
vith  |s|  ^1,  every  partition  Sl  and  every  set  of  distinct 

values  with  0 < |kj^1  < 1, 

(3«10)  det[ B( a> . . * ,B(z> Kjj»n^) 3 ^ 0. 

Proof.  Take  an  arbitrary  with  |zq|  >1  and  let 
1 < a < k,  be  the  distinct  roots  of  (3.1),  each  with  multiplicity  »q(*q). 
Since  (“a(*o)^^l  * partition  of  r,  then  (3.10)  holds  for  our  z^, 

"0^*0^  Lemma  3.1,  z^  is  not  an  eigenvalue  of  G, 

and  by  Kreiss  Main  Theorem  in  [4],  stability  follows. 

Since  the  determinant  in  (3.10)  depends  entirely  on  the  boundary  functions 
(3.7),  Theorem  3.2  is  scheme-independent.  Thus  in  applying  the  theorem, one 
avoids  the  inherent  difficulty  of  solving  the  characteristic  equation  (3.**). 

4.  Translatory  Boundary  Conditions 

In  this  section  we  return  to  consider  the  approximation  (2.4) (2.6)  where 
the  boundary  conditions  are  translatory. 

The  boundary-functions  associated  with  (2.6)  are 
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b (B)K)  • 2 (1  ■ Ofaxff'l* 

^ J-0  J ^ 

niua  writing  b(z,K)  ■ bQ(zfK)  we  bzw 

b^(z,K)  • b(Z(K)y  V ■ Oy***}!*'!* 

n 

Consequently,  If  1*  * partition  of  r,  then  the  rxr  matrix 

(4.1)  (B(z,Kj^,m^),...,B(z,..^,mjj)l 


In  (3.10)  is  given  by  the  rxmj^  blocks 


b(z,K) 

b(z,K) 

b(s,K) 

B(z,ic^,mj^)  « 

K b(z,K) 

• 

• 

d 

K b(z,K) 

• 

• 

^ 

K b(z,K) 

• 

• 

m 

^(z,k) 

k*‘"^(z,k) 

a*  m 

K**  ^(z,k) 

The  fact  that  (4.1)  is  determined  now  by  the  single  boundary- function 
b(z,K),  Implies  the  following  significant  simplification  of  Theorem  3.2. 


THEOREM  4.1.  Approximation  (2.4) (2.6)  is  stable  if  for  • 
every  z with  |z|  > 1 and  k with  0 < |k|  < 1,  we  have 

(4.2)  b(z,K)  . 2 (zc$^)  - c(0))kJ  4 0. 

J-0  ^ ^ 


Proof.  Take  an  arbitrary  z with  |z|  >1,  a partition 
of  r,  and  distinct  values  k^,  ^ ^ with  0 < |k^|  <1.  In 

order  to  prove  stability,  it  suffices,  by  Theorem  3.2,  to  verify  (3. 10). 
For  this  purpose,  let 


r-1 

U)»0  ■ 


.1"^ 

(i'iJ  '>(*»'^)I 


- 0, 
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be  a vanlshlxig  linear  combination  of  the  rows  of  (4.1).  The  vector 
relation  in  (4.3)  consists  of  r scalar  equations, 


[ 


^ r„uv/.^,)]  j 


T y b(z, 

«o  V 


Ui?0 

vhlch  we  write  as 


K>Kj 


» 0,  l<i<N,  0<J<mj^-l, 


(4.1.)  .5iL,„)r  Y «“1|  . 

I L n=o  ^ J ^ tt=Kj^ 


l<i<N,  0<J<mj^-l. 


Since  0 < |k^|  <1,  then  by  hyi)othesis,  the  expression  in  the  left 
brackets  of  (4.4)  satisfies 

^^(**'')^k=Kj  / 0,  1 < i < N. 

Thus,  expanding  (4.4)  by  Leibniz'  rule  and  using  induction  on  J > 0, 
we  find  that  the  sum  in  (4.4)  have  vanishing  derivatives  at  k = k i.e.. 


r u 1 

1 

1 ” 1 ' “> 

Lu=0  J K=Kj^ 

l<i<N,  0<J<m^-l. 

We  conclude  that  the  polynomial 


P(»c)  fi  E y y 
u=0  ^ 


which  is  of  degree  r-1  at  most,  has  r roots;  k^,  1 < i < N,  each  with 
Bultlpllclty  mj^.  Hence,  P(k)  e 0 and  the  coefficients  must  vanish. 
By  (4.3),  therefore,  the  rows  of  (4.1)  are  linearly  Independent,  so  (3.10) 
holds,  and  stability  follows. 

Before  turning  to  the  proof  of  the  Main  Theorem  (Theorem  2.1)  we  note 
that  Theorem  4.1  applies  also  to  single-leveled  boundary  conditions.  For 
example,  we  can  Immediately  obtain  the  following  result  of  Kreiss  [3]  (see 
also  [1]):  Let  the  boundary  conditions  be  determined  by  extrapolation  of 
arbitrary  degree  s - 1,  i.e., 


IS 


(k,3) 

The  associated  boundary*function  is 

b(K)--i;  (®)(-1)^kJ  = -(1‘-k)®. 

J=0  ^ 

Thus  b(K)  / 0 for  0 < |k|  <1,  emd  by  theorem  V.l,  (2. 4) (4.^)  is  stable. 

Proof  of  Theorem  2.1.  The  amplification  factor  T(|)  of  the  boundary- 
scheme  (2.7)  is  given  by 

f(0  = a =0,1. 

By  (2.8)  we  have  f^^^(5)  / 0,  thus  $(|)  is  well  defined.  Since  the 
boundary- scheme  is  stable  then  |f(|)|  <1,  so  finally 

(4.6)  l*<®>(Ol  < /O  V6. 

The  boundary-function  associated  with  (2.6)  satisfies 

b(z,K)  = U (*c(^^  - c^®^)k^  c zT^^\k)  - T^®^(k). 
d»o  ^ ^ 

So,  for  |z|  > 1,  we  use  (4.6)  to  find  that 

£ Ul  • |T<^>(!)I  - >0. 

That  is,  the  equation 

(4.7)  1>(z*k)  = 0 with  |z|  >1, 

has  no  roots  k with  |k|  * 1.  Since  the  roots  of  (4.7)  are  continuous 
fhnctions  of  z,  we  conclude  that  the  number  of  k with  |k|  < 1 is 
fixed  for  |z|  >1,  and  can  be  detezmlned  by  considering  large  values  of 
|s|.  Writing  (4.7)  in  the  form 


T^^^(k)  - z“V°\k)  . 0, 


ve  let  |z|  -•<»  and  use  (2.8)  to  find  that  (4.7)  haa  no  roots  in  .tait 
disc.  In  other  words,  if  |z|  > 1 and  b(z,K)  « 0,  than  k Mut  satisfy 
|k|  > 1.  By  continuity  therefore,  if  |z|  > 1 and  b(s,K)  > 0,  then 
|k|  > 1}  i.e., 

b(z,K)  0,  N/|z|  >.  1,  |k|  < 1» 

This  ioQ>lies  (4.2),  and  Theorem  4.1  completes  the  proof. 
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